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ABSTRACT 

We use cosmological SPH simulations to study the kinematic signatures of cool gas accretion onto a pair of well- 
resolved galaxy halos. Cold-flow streams and gas-rich mergers produce a circum-galactic component of cool gas that 
generally orbits with high angular momentum about the galaxy halo before falling in to build the disk. This signature 
of cosmological accretion should be observable using background-object absorption line studies as features that are 
offset from the galaxys systemic velocity by ~ 100 km/s. Accreted gas typically co-rotates with the central disk in the 
form of a warped, extended cold flow disk, such that the observed velocity offset is in the same direction as galaxy 
rotation, appearing in sight lines that avoid the galactic poles. This prediction provides a means to obscrvationally 
distinguish accreted gas from outflow gas: the accreted gas will show large one-sided velocity offsets in absorption line 
studies while radial/bi-conical outflows will not (except possibly in special polar projections). This rotation signature 
has already been seen in studies of intermediate redshift galaxy-absorber pairs; we suggest that these observations may 
be among the first to provide indirect observational evidence for cold accretion onto galactic halos. Cold mode halo 
gas typically has ^3 — 5 times more specific angular momentum than the dark matter. The associated cold mode disk 
configurations are likely related to extended HI/XUV disks seen around galaxies in the local universe. The fraction of 
galaxies with extended cold flow disks and associated offset absorption-line gas should decrease around bright galaxies 
at low redshift, as cold mode accretion dies out. 

Subject headings: cosmology: theory — galaxies: formation — galaxies: halos — methods: TV-body 
simulations — methods: SPH 



1. INTRODUCTION 

Arguably the most substantive change in our theo- 
retical picture of galaxy formation in the last decade 
has involved the deposition of gas into galaxies. Cos- 
mological simulations set within the LCDM frame- 
work show that the majority of galactic baryons 
should be accreted in a "cold mode", which docs 
not shock hea t in the ha l o pri o r to building the 
galaxy itself (IKeres et al.l |2005 | ; iDekel fe Birnboim 
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20091: IDekel et alJ 120091: IKeres et al 
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20111: IStewart et al 



Ivan de Voort et al.1 1201 lF iFaucher-Giguere et al 
201 If ). This is quite different from the canonical pic- 
ture where gas first obtains the halo virial temper- 
ature and spin before fallin g in to fo r m an angular 
momentum supported disk dSilld 1 19771 iWhite fe ReeT 



19781 iBarnes fc Efstathioul 119871 : IWhite fe Frenkl 119911 
Mailer fe Bullock! 120041 ) . Instead, short cooling times 



prevent the development of a stable shock near the virial 
radius in halos that ar e less massive than a cr itical 
threshold (|Binnevl \WJ% iBirnboim fe Dekell [2003]) and 
cool gas gets deposited directly into galactic halos. 
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Previous simulation studies have noted that this gas 
seems to be rotating with large amounts of angu- 
lar momentum, which may help explain what hap- 
pens to t he significant poten t ial energy of the infalling 
cool gas (IKeres et all I2009bt IKeres fe Hernquistl 120091: 
lAgertz et all 120091 : iBrook et al.1 120111 ). Here, we study 
the cold mode gas and its detectability in detail, demon- 
strating that this accreted cool gas orbits about the halo 
before falling in to build the central disk, delivering not 
just fuel for star formation but also angular momentum 
to shape the outer galaxy. This angular momentum en- 
ables an observable prediction: cosmological accretion 
should be discernible in halo absorption-line studies as 
cool components that are offset from the galaxy's sys- 
temic velocity by ~ 50 to 200 km s -1 , co-rotating with 
the central disk in most cases. 

Though the cold mode accretion prediction seems to 
be fairly robust and well understood analytically, there 
remains no smoking-gun observational evidence that cold 
mode accretion is actually occurring in nature. It is im- 
portant to look for such a signature because it is always 
possible that additional physics not currently included 
in simulations changes the story. Such a situation is 
not out of the question, given that the same simulations 
that are used to predict and characterize cold mode be- 
havior also over-predi ct the baryonic ma ss function of 
galaxies significantly (|Keres et al.l l2009at ). In order to 
reconcile the relatively high level of baryonic deposition 
into galaxies with a correspondingly small fraction of 
baryons locked up i n galaxies today (|Guo et al.1 120101 : 
iBehroozi et alJ 120101 ) , one must appeal to outflows of 
some kind. Indeed, observations of cool halo gas around 
galaxies (typically via metal line absorption systems) 
have em phasized the presen c e of gas outflows from galax- 
ies (e.g., ISteidel et alJll996l: IShaolev et al.ll2003l: IMartinl 



127)051 : 1 Weiner et al.l 2009: Steid el et al.H2010l : IRubin et al 
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The clear observational evidence for outflow (rather 
than inflow) is at least understandable at high rcdshift; 
gas accretion at z > 2 is expected to flow along particu- 
larly dense filaments that can penetrate into the center of 
halos, even above the aforementioned critical halo mass. 
Dense, narrow streams of this kind may result in glob- 
ally s mall covering fractions, making them difficult to ob- 
serve (iFaucher-Giguere k Kerell20lll: iKhnm et al.ll20lH: 
Fuma.gal het al.ll2011[ ). Galaxies at z ~ 2 are also at the 
peak of cosmological star formation (|Hopkins k Beacoml 
120061 ): at these epochs, star-formation driven feedback 
processes are even more likely to obscure any indicators 
of gas accretion (though feedback processes will invari- 
ably obscure gas accretion signatures at all epochs, to 
degrees that are largely unknown). Nevertheless, even 
at moderate redshifts [z < 1.5) there are no definitive 
observational indications of significant cosmological gas 
accretion into galaxies in the form of cold mode accre- 
tion. 

In a companion paper (|Stewart et al.|[20Tl|) we showed 
that while the covering fraction of cool accreted material 
remains relatively high in low-mass halos at moderate 
redshifts (~ 30 - 50% within R < 50 kpc at z ~ 1), the 
covering fraction of such gas should drop sharply in halos 
more massive than the critical threshold for sustaining 
virial shocks (M v j r ~ 10 12 M Q ). In principle, one would 
like to examine quasar-galaxy (or galaxy-galaxy) pairs 
over a range of galaxy masses in order to witness the 
end of the cold flow regime galaxy-by-galaxy. But in 
order to do this, one must differentiate absorption signals 
from accreted gas and absorption from gas that has been 
deposited there from outflows. In this paper we show 
that the kinematics of absorbing gas can provide a means 
to distinguish its class: halo gas from cold mode accretion 
generally rotates about the galaxy. 

It should not be surprising that cold-mode gas tends to 
rotate with the disk. This is the primary mo de by which 
disk galaxies are built in simulations (e.g., iDekel et al.l 
12009( 1. Cold flows and gas- rich mergers not only supply 
mass to galaxies but also their angular momentum. In 
fact, we find that cool accreted gas tends to have higher 
and more coherent specific angular momentum than the 
dark matter in the halo and this is one property that 
enables an observable kinematic signature. 

In the next section we describe our two simulations, 
each of which tracks the evolution of a Milky Way-size 
galaxy down to z = 0, and we detail our method for 
creating mock absorption sight lines through the simu- 
lations. We briefly explore the general properties of the 
cool gaseous halos of our simulated galaxies in Sj3] In Sj4] 
wc analyze one of the simulated galaxies at a single point 
in time, detailing a unique observational signature asso- 
ciated with orbiting halo gas, in the form of an extended 
disk of cool gas that co-rotates with the galaxy. We in- 
vestigate the statistical significance of this result Sj5l an- 
alyzing both simulations over a broad range in redshift. 
We compare this co-rotation signature to observations of 
galaxy-absorber pairs, and discuss the implications of our 
results in S|6] We summarize our findings and conclude 
in g3 

2. SIMULATIONS AND ANALYSIS 



2.1. Simulations 

Each of our two simulations utilize a separate set of 
cosmological initial conditions, tracking the evolution of 
a Milky Way size dark matter halo until z = 0. Be- 
cause we primarily investigate the single most massive 
galaxy in each simulation, we refer to these simulations 
throughout this paper as "Halo 1" and "Halo 2." Halo 1 
has an active merger history until z ~ 1.5, but is subse- 
quently quiescent (a dark matter only simulation of the 
same initial conditions was performed at very high reso- 
lution in the Via Lactea II simulation of Dicman d et al.l 
(|2008ll 5 l ). In contrast, Halo 2 experiences a relatively qui- 
escent early history, but has a major merger at z ~ 0.8, 
and another large merger at z ~ 0.2. For both simu- 
lations, the most massive galaxy has a dark halo virial 
mass of M vir = 1.4 x 1O 12 M at z = 0. 

We use "zoom-in" multi-mass particle grid ini- 
tial conditions gen erated with the GRAFIC-2 package 
(|Bertschin gcr 2001). A periodic box of 40 co-moving Mpc 
is used in each simulation to properly account for large- 
scale tidal torques. Each simulation contains a high res- 
olution region, limited to a ~ 6 co-moving Mpc cube; in 
this region, the masses of the simulation particles in the 
initial conditions are (mdark, «^ g as) = (17, 3.7) x 1O 5 M Q 

We adopt the best-fit cosmologi cal parameters of th e 
WMAP three-year data release (jSpergel et al.l 12007ft : 
fl M = 0.238, fl A = 0.762, H = 73km s" 1 Mpc" 1 , 
n s = 0.951, and cr 8 = 0.74 

We use the code GASOLINE (|Wadslev et. a]ll200l . 
which is a smoothed particle hydrodynamics (SPH) ex- 
tension of the pure T V-Bod y gravity code PKDGRAV 
developed by iStadell (|200lD . We use a gravitational 
force softening of 332 pc, which evolves co-movingly un- 
til z = 9 and remains fixed from z = 9 to the present. 
The SPH smoothing length adapts to always enclose the 
32 nearest gas particles and has the minimum set to 0.05 
times the force softening length. Renderings of the two 
simulations are shown in the upper (Halo 1) and lower 
(Halo 2) panels of Figure[TJ The panels from left to right 
show progressively zoomed renderings of the gas density 
in each simulation at z = 2. 

The code assumes a unifor m UV backgroun d from 
QSO, implemented following lHaardt k Madaul (|1996l ) 
and F. Haardt (2002, private communication). It in- 
cludes star form ation and Compto n and radiative cooling 
as described in iKatz et all (|1996l) . The code calculates 
the abundance of neutral hydrogen by assuming an op- 
tically thin ideal gas of primordial composition and in 
ionization equilibrium with the UV-background, treat- 
ing collisional ionization, photoionization and recombi- 
nation pro c esses. Rates for collisional ionizatio n from 
lAbel et all (119971). radiative recom bination from iBlackl 
(|1981[ ) and lVerner k Ferlandl (|1996T > have been used. 

In this paper, we focus on cool gas accretion and the 
possibility of detecting it as quasar absorption systems, 
which primarily utilize metal lines (Mgll, CIV, OVI, 
etc.). To correctly produce metal lines in a simulation 
of this kind requires radiative transfer, metal diffusion 

8 The merger history and other data about this halo can be 
found at http://www.ucolick.org/~diemand/vl/ 

9 For c omparison, recent simula t ions b y [Kcrcs & Hc rnauisH 
(2009) and Fauchcr-Gigucrc k, KcreJ ( f20TTTl have particle masses 
of 9 X 10 4 M Q and 4 X 10 4 M Q , respectively. 




Fig. 1. — A visualization of projected gas density in each of our simulations at z = 2, with the color variations on a log scale with red, 
green, and blue corresponding to gas column densities of ~ 0.03, 0.4, 2.0 Mq/pc 2 , respectively. The top panels show Halo 1 while the 
bottom panels show Halo 2. From left to right, panel widths correspond to 5 Mpc, 1 Mpc, and 200 kpc in co-moving coordinates. Note 
that both of our main galaxies reside within an extended filamentary structure at this redshift. 



and modeling of local ionizing sources, which are not in- 
cluded here. Since we are only focusing on the qualitative 
behavior of halo gas, we will instead give results in terms 
of HI column density, calculated in the optically thin 
limit, without local sources. For column densities below 
the Lyman limit (2 x 10 17 cm~ 2 ) our computed column 
densities should be fairly robust. For higher column den- 
sities, we expect a full treatment to lead to quantitative, 
but not qualitative differences. (See £13.11 for a quanti- 
tative comparison of our covering fractions with th e re- 
cent simulations bvlFaucher-Giguere fc Kc rcs (|2011[) and 
iFumaealli eTail (|20il[ )" that do contain radiative transfer 
treatments.) 

The "blastwave" supernova feedback model imple- 
mented in our simulations creates turbulent motions in 
nearby gas particles that kee ps them from coolin g and 
forming stars, as described in lStinson et. all ([2006D . The 
only free parameters in the star formation and feedback 
model (minimum density threshold, p mnl = 0.1 cm~ 3 ; 
star formation efficiency factor, c* = 0.05; and fraction of 
supernova energy that co uples with the ISM, eSN = 0.1) 
have been motivated by IGovernato et al.l (|2007f ) in or- 
der to produce galaxies with a realistic star formation 
rates, disk thicknesses, gas turbulence, and Schmidt law 
over a range in dynamic masses. The feedback model 
we use here is similar to those used in recent simula- 
tions that have shown great success in producing real- 



istic disk- type galaxies (|Governato et al.l 12009ft. match - 
ing the mass-metallicitv relation ([Brooks et al.l 120071) , 
and matc hing the abundance of Damped Lya systems 
at z = 3 (IPontzen et all 120081 ). We refer the reader to 
IGovernato et al.l (|2007tT for a more detailed description 



of the simulation code. 

With the resolution implemented here, this feedback 
model results in minimal winds (~ 100 km/s) that mostly 
affect hot gas, a nd are more prom inent in halos with 
M vir < lO n M (jShen et al.|[20Tot) . Although outflow- 
ing gas may be more prominent in higher resolution 
simulations with higher density thresholds and modified 
star formation and sup ernova efficiency parameters (e.g., 
IGovernato et al.ll2010f) , we emphasize that the halos dis- 
cussed here do not have outflows of cool gas that would 
be detectable through absorption. 

We define the viri al radius for our halos by 
iBrvan fc Norman] (|1998|) , noting that this is a fairly typ- 
ical definition used in constructing merger histories from 
jV-body simulation s of dark matter substructure (e.g. 
iStewart et al.l 120081 ). While this paper primarily focuses 
on mock observations of absorption lines through our 
simulated galaxy halos, we analyze the formation and 
evolution of the galaxies themselves in greater detail 
in an upcoming paper. Both galaxies are disks before 
they reach the transition mass for hot mode accretion 
(M V j r ~ 10 12 Mq). After this transition, both galaxies 
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Fig. 2. — Example absorption sight lines taken from Halo 2 at 
z=1.4, showing HI column density per unit velocity versus line of 
sight velocity. In each panel, the galaxy systemic velocity is 
km/s, indicated by the vertical dashed line. Left: sight lin es th at 
show orbiting gas, with velocity offsets from systemic (see ^12.21 for 
definitions). The top and middle plots are co- rotating, due to their 
spatial positions with respect to the galaxy, while the bottom panel 
is anti-rotating (i.e., the central galaxy's rotation curve would be 
positive on this side of the galaxy). The impact parameters of 
these particular sight lines are (from top to bottom): 50, 20, 60 
co-moving kpc. Right: sight lines that are not orbiting the galaxy, 
but show velocities near (or on both sides of) systemic. The impact 
parameters of these sight lines are (top to bottom): 6, 20, 40 co- 
moving kpc. 

lose their reservoirs of cool high angular momentum halo 
gas. Halo 1 grows a massive bulge as hot gas cools from 
the halo, seemingly due to secular instabilities (though it 
is possible that numerical resolution issues may give rise 
to artificial bulge-creation in the central regions, as well). 
Halo 2 develops into a more spheroidal system due to a 
major merger just prior to reaching the transition mass. 
We note that previous works that have produced disk 
galaxies at z = using a similar code as the one used 
here typically involve halos less massive than those w e 
study here (e.g., 7 x 1O U M in lGovernato et al.ll2009h . 
which may be the cause of these morphological differ- 
ences, since more massive galaxies are known to host a 
small er fraction of disk galaxies (e.g. IWeinmann et al.l 
[2009ft . 

2.2. Mock Absorption Sight Lines 



In order to study the kinematics of cool halo gas in a 
manner useful for comparison to observations, we con- 
struct regularly spaced lines of sight around our galaxies 
for a variety of viewing angles and for different times. 
For each sight line, we use the analysis software TIPSY's 
"absorb" command to compute the column density of 
neutral hydrogen per unit velocity, as a function of ve- 
locity along the line of sight, mimicking the practice of 
searching for absorption gas along a quasar line of sight. 
The column density is computed by first dividing each 
sight line into spatial bins of order the size of the spa- 
tial resolution of the simulation. For each spatial bin, 
each gas particle with a smoothing length that intersects 
the line of sight contributes to the total mass density of 
the bin. For each of these particles, the contribution of 
the spline kernel to the column density of each bin that 
overlaps the kernel is integrated along the line of sight. 
The velocity of each bin is the mass-weighted velocity 
of particles contributing to the total column density of 
neutral hydrogen. Each of these bins are broken into a 
number of sub-bins, with gas velocity and temperature 
linearly extrapolated from adjacent bins. Velocity pro- 
files for each sub-bin is based on thermal broadening and 
dampening wings, with each sub-bin's profile integrated 
over the velocity bin to determine its contribution. 

We limit our examination to absorption lines with 
iVm > 10 16 cm -2 , which allows us to focus on the fairly 
low column densities that our simulations are most capa- 
ble of reproducing]- This range also facilitates our later 
comparison with weak Mgll absorbing gas, as it is ap- 
proximately the mini mum HI column densi t y associated 
with Mg II detections ([Churchill et al.ll200Ct iRigbv et al.l 
12001 . 

Figure [2] shows six example sight lines from our simu- 
lations (for Halo 2 at z = 1.4) that meet our iVm > 10 16 
cm -2 criterion. In all panels, we plot the column den- 
sity of neutral hydrogen per velocity bin as a function 
of velocity, with the galaxy systemic velocity set to 0, 
represented by the vertical dashed line. In order to com- 
pare the kinematics of the absorbing gas to the central 
galaxy, we classify each line of sight as either orbiting or 
not orbiting. 

We define a sight line to be orbiting if > 90% of the 
absorbing gas lies entirely to one side of the systemic ve- 
locity of the galaxy, where the galaxy's systemic velocity 
is the Hi-weighted projected velocity through the center 
of the galaxy. Examples of sight lines classified as or- 
biting are shown by the left three panels in Figure [2] If 
this criterion is not met, then we classify the kinematic 
signature as not orbiting. Note that this is a fairly con- 
servative classification, in the sense that we classify even 
ambiguous cases as not orbiting. In the top right, the 
classification of not orbiting is due to the presence of two 
distinct components of gas that straddle the systemic ve- 
locity (Interestingly, this configuration is reminiscent of 
something that might be tagged as "outflow" in an em- 
pirical study, though in our case, it is simply the chance 
intersection of two clumps of accreted gas that happen 

10 We again emphasize that comparing directly to observed 
metal lines would require radiative transfer, metal diffusion and 
modeling of local ionizing sources, which are not included here. 
This is why we instead focus on the properties of HI gas at column 
densities below the Lyman limit (2 X 10 17 cm -2 ). 
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to be orbiting in different directions.) Such a configura- 
tion for is rare, however; more common examples of gas 
that is not orbiting are the other right panels, each of 
which has a single distribution of halo gas near the sys- 
temic velocity. In contrast, the left middle panel shows 
some absorption on both sides of systemic, but because 
> 90% of the gas lies to a single side, this sight line is still 
considered orbiting. This configuration is also rare. Sin- 
gle peaked, relatively narrow absorption lines as shown 
in the two bottom panels are the most common profiles, 
with > 95% of sight lines having a full width half maxi- 
mum velocity spreads < 110 km/sP"! 

We note that in detail orbiting gas may not necessarily 
show a velocity offset from systemic along an arbitrary 
projected line of sight. Also, gas that does show a veloc- 
ity offset from systemic will have unknown radial veloci- 
ties that could be substantially greater than the velocity 
offset detected. Despite these shortcomings, we choose 
to characterize these velocity-offset sight lines as "orbit- 
ing" because a full 3-dimensional analysis of our galaxies 
demonstrates that cool halo gas detected by a velocity 
offset from systemic is typically orbiting the galaxy in a 
relatively coherent fashion (see 

If a sight line is found to be orbiting, we further clas- 
sify it as either co-rotating or anti-rotating by comparing 
the direction of the velocity offset with the line-of-sight 
rotation of the galactic disk along the same projection 
(sight lines to absorbing gas at near right angles to the 
galactic disk arc not included in this classification, as pro- 
jection effects make these sight lines ambiguous in terms 
of co-rotation or anti- rotation). As an example, in the 
left panels of Figure [5J the top and middle panels are co- 
rotating, due to their projected position with respect to 
the galaxy, while the bottom left panel is anti-rotating. 

3. HALO GAS PROPERTIES 

Before using our simulated galaxies to probe kinematic 
signatures of cool halo gas, it is important to know the 
general properties of our galaxies and their gaseous ha- 
los. Because our primary focus is cool halo gas seen in 
absorption, we choose to first investigate the radial ex- 
tent of cool gas, to define a useful outer radius within 
which to focus further analysis. 

3.1. Covering Fraction Versus Radius 

Some past studies of metal line absorption systems as- 
sociated with galaxies suggest the presence of a critical 
radius for cool halo gas. In this model, absorption gas ex- 
tends in a roughly spherical configuration out to a given 
radius, and gas is unable to cool beyond this radius (e.g 



ISteidellH995TlTinker fc Chenll200a IChen fc Tinkei1l2008h . 
Other studies suggest that the covering fraction of de- 
tectable a bsorption systems should decreases as a po wer 
law in R (jMartin fc Bouchell2009l : iSteidel et al1[2010h . 

In order to quantify the radial extent of cool halo gas, 
we define the covering fraction of accreted neutral hydro- 
gen as a function of radius, f c (<R), as the total fraction 
of sight lines (within a projected radius R from the cen- 



ter of the galaxy) for which Nm > 10 cm (see 
for details on absorption sight line construction). Note 

11 Because our feedback model does not produce cool gas out- 
flows, our absorption sight lines are typically within 500 km/s of 
systemic, and never offset by more than 1,000 km/s. 
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Fig. 3. — Example covering fraction versus projected radius pro- 
files for Halo 2 at various epochs. The radial profile of f c (<R) is 
similar across a broad redshift range, so long as the galaxy is not 
experiencing a major merger. The clashed green line shows a power 
law fit given by EquationTTl 

that we include the galaxy disk in this analysis, so that 
f c ~ 1 for small values of R. At large radii, this inclusion 
of the galaxy to the total covering fraction is relatively 
minor (< 10%, 4%, 1% of sight lines for R < 30,50,100 
co-moving kpc). 

Figure [3] shows five example profiles (for Halo 2 at 
z ~ 3, 2, 2.5, 1.3, 1) of the covering fraction as a function 
of radius, from R < 10 to R < 100 co-moving kpc, aver- 
aged over three orthogonal orientations of the galaxy, for 
Njzi > 10 16 cm~ 2 . Our simulated galaxies do not show 
a distinct truncation radius that cleanly separates a cold 
gas regime, but instead we find a smooth progression 
from higher covering fractions at small radii, to smaller 
covering fractions at larger radii (though this profile is 
less smooth during a gas-rich merger). As shown in Fig- 
ure [3j our f c vs. R profiles are quite similar over a wide 
redshift range. Each curve is well fit by a power law in 
R, within the regime 10 < R < 100, with R in co-moving 
kpc: 

Rq 



.U<R) = 



(1) 



where Rq represents an inner radius within which the 
covering fraction is ~ 100%. The dashed green line in 
Figure [3] is given by Rq = 10 co-moving kpc and j3 = 0.7. 
We note that our result s are in qualitative agreement 
with ISteidel et al.l (|2010D . who found the covering frac- 
tion around LBG galaxies at z ~ 3 to decline as a power 
law in R, with a slope of 0.2 — 0.6, very similar to our 
value of 0.7. In addition, our results are also qualita- 
tively simil a r to t he lower redshift sample (z < 0.5) of 
IChen et al.l (|2010D . who found that the covering fraction 
is near unity at impact parameters R < 20 kpc, and de- 
clines roughly linearly in log R to < 5% covering fractions 
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at R ~ 100 kpc. 

In detail, the profile we present above also depends on 
the minimum value of -ZVhi, with steeper radial profiles 
for higher column density gas (similarly, higher equiv- 
alent with absorption system s also show steeper radial 
profiles, see lChen et al.lf^OlOD . Accounting for both sim- 
ulated galaxies, we find more generally that the covering 
fraction as a function of radius and minimum column 
density (from Ahi > 10 13 cm" 2 to Am > 10 18 cm" 2 ) is 
well-fit to the following functional forrrFl: 

/ R \ _,3(ArHl) 
f e (N m , <R)= , (2) 

P(N m ) = Pie + 0.1 log[A^Hi/10 16 cm- 2 ] (3) 

with i? ~ 10 co-moving kpc and /3i6 ~ 0.6—0.8. Because 
the cool gaseous halos of our galaxies die out after cold 
mode accretion ends (z < 1.3 and z < 0.8 for Halo 1 
and Halo 2, respectively) we caution that this derived fit 
for / c (<i?) is based on a narrow halo mass and redshift 
range: M vir ~ 10 11 - 10 12 M o and 3 < z < 1. 

We also caution that our simulations do not have treat- 
ment of radiative transfer (see §2.2j) . which may affect 
the detailed values of AThi presented here. In order to 
give some estimate of how this may impact our result- 
ing covering fractions, we compare the covering fraction 
of our galaxies at z = 2 with rece nt hi gh-z simulations 
of IFaucher-Giguere fc K crcs l|2011[) and IFumagalli et aTl 
(|2011[ ). both of which include post-processing radiative 
transfer prescriptions for galaxies of similar mass to 
those analyzed here. At z ~ 2, these works both find 
that / c (< 2i? vir ) - 4% and f c (< i? vir ) - 10% for 
Ahi > 10 17 ' 3 cm" 2 , while our simulations (utilizing the 
above best-fit equation, and i? v i r ~ 250 co-moving kpc 
at z = 2) have covering fractions of ~ 4% and ~ 7%, 
respectively, for the same column densities. 

This indicates that the lack of radiative transfer may 
lower the available HI content in the halos of our galaxies, 
suppressing the covering fraction by a factor of ~ 30% at 
a fixed column density threshold. (This, in turn, suggests 
that our reported column densities may be slightly higher 
than they would be with inclusion of radiative transfer) . 
Because the rest of this paper primarily discusses the 
kinematic properties of cool halo gas in our galaxies, only 
utilizing the integrated column densities of sight lines to 
determine if they are potentially observable or not, we do 
not believe the lack of radiative transfer will substantially 
alter our findings. 

It is worth noting that most of the accreted halo gas 
resides within R < 100 co- moving kpc of our simulated 
galaxies, with the only notable exception being major, 
gas-rich galaxy mergers that have not yet fallen within 
this radius. (Mergers of this kind are more likely at high 
redshift, when mergers are more common, and galax- 
ies are more gas ric h than at lower redshift; see e.g., 
I Stewart et al.l l2009alTbh . Since the infall time of m ajor 
mergers is short (e.g., IBovlan-Kolchin et al.l 120081 ) . we 
will focus our analysis on R < 100 co-moving kpc for the 
remainder of this paper. For our simulated halos, this 

12 We note that under the presumption of a power law form for 
f c (<R), the differential covering fraction of gas precisely at the 
projected radius R must follow the same power law slope, with a 
normalization that is lower by a factor of (1 — 13/2). 



choice of radius corresponds to ~ i? v ; r /2 at z = 3 and 
~ -Rvh /3 at z = 0. 

3.2. Mass and Angular Momentum Evolution 

The top panels in Figure 0] show the mass buildup of 
each galaxy over cosmic time, with the left and right 
panels showing Halo 1 and Halo 2, respectively. In each 
panel, the solid black line shows the total virial mass of 
each halo (including dark matter and all baryons within 
R < R vn ), while the other curves represent the baryonic 
masses within R < 100 co-moving kpc; the dashed green, 
dot-dashed red, dotted blue and dotted magenta lines 
show the stellar mass, hot gas mass (T > 10 5 K), cool 
gas mass (T < 10 5 K), and HI gas mass, respectively. 

Once each galaxy crosses the critical halo mass for 
forming shocks, M s h, accreted gas is shock heated 
strongly enough that very little cool gas reaches the in- 
ner halos of our galaxies. This point is indicated by ver- 
tical lines in each panel of Figure [4] We see that after 
this critical transition, the gas mass in the halo becomes 
dominated by the hot component (red dot-dashed line) 
and the cool component drops in comparison (blue dot- 
ted line). At the sa me time, the coverin g fraction of cool 
gas drops severely ([Stewart et al.ll20Tll) and the galaxies 
begin to evolve into more bulge-dominated morphologies, 
rather than disks (see £16.31 for more on this transition). 

Given that the velocity distribution of halo gas is the 
key observational signature we are exploring, it is use- 
ful to characterize the kinematics of orbiting halo gas in 
terms of its angular momentum. A useful way to do thi s 
is via a dissipationless spin parameter (|Peeblesl I1969T) . 
Specifically, we compare the specific angular momentum 
j of the dark matter, hot gas, and cool gas as a function 
of time in eac h halo using the spin parameter defined by 
IBullock et all (|200l 



Here j x is the specific angular momentum of component 
x (= dark matter, hot gas, cool gas, and all gas) within a 
sphere of radius R and V is the circular velocity measured 
at radius R. 

The bottom panels of Figure HI show the spin param- 
eter for the dark matter within R = 100 co-moving 
kpc as a solid black line. The spin parameter evolu- 
tion of the gaseous components within the same radius 
is represented by the black dashed (all gas), blue dot- 
ted (cool gas), and red dot-dashed lines (hot gas). Note 
that the spin parameter in dark matter evolves some- 
wha t stochastically, spinning up in response to merg- 
ers (jVitvitska et al.l 120021: IMaller et ail I2002D but then 
settling down as the halos re-virialize and the angula r 
momentum is redistributed (jD'Onghia fc Navarroll200"7[) . 
The redistribution does not happen as readily in gas, 
which tends to cohere and spin up more strongly. The 
average spin for the dark matter bounces around Adm ~ 
0.04 over time, at a value that is quite similar to the ex- 
pectations of many past studies from dissipationless sim- 
ulations (recently iBett et al l [20071 IMaccio et all 120071: 
iBerta et al.ll2008fc iBett et al.ll2010f T On the other hand, 
the gaseous halo components have much more angular 
momentum, demonstrating a more sustained response to 
mergers and inflow. The cool gas has the highest spin 
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Fig. 4. — Baryonic masses and spin parameters (R < 100 co-moving kpc of each simulated galaxy) versus time. The left panels show 
Halo 1, while the right panels show Halo 2. Top: the total mass within the virial radius is shown by the solid black line, while the stellar 
mass, hot gas mass (T > 10 5 K), cool gas mass (T < 10 5 K), and HI gas mass within R < 100 co- moving kpc are given by the dashed 
green, dot-dashed red, dotted blue, and dotted magenta lines, respectively. Bottom: the spin parameter of the dark matter and gas within 
R < 100 co- moving kpc (see Equation |4j . The dark matter spin parameter is given by a solid black line, while the spin parameters of 
cool gas, hot gas, and all gas are given by the blue dotted, red dot-dashed, and black dashed lines, respectively. Note that after rapid 
accumulation of angular momentum from mergers, the dark matter spin parameter quickly settles to pre-merger values, while the gas spin 
parameters stay higher considerably longer. The vertical line in each panel shows where t he galaxy crosses the transition from cold mode 
to hot mode accretion, M v i r = Af s h ~ 10 12 Mq, as motivated by Dckcl & Birnboim (2006) for shocks at a moderately large fraction of the 
virial radius. After this threshold, little cool gas reaches the galaxy from cold flow accretion. 
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Fig. 5. — Line of sight velocity of cool halo gas in Halo 2's main 
galaxy at z = 1.4. The image width is 200 co- moving kpc, and the 
galaxy is viewed near edge-on at an inclination of 75. The black 
contours enclose regions with JVhi > 10 16 cm -2 , and the blue-red 
shading shows the average (mass- weighted) velocity of HI along the 
line of sight, from —300 to 300 km/s. The galaxy (inner region) 
tends to rotate in the same direction as cool halo gas, even out to 
100 co- moving kpc). 



of all, some ~ 3 — 5 times higher than that of the dark 
matter with A coo igas ~ 0.1 — 0.2. Wc emphasize that 
this result is particularly important for theoretical stud- 
ies that try to understand angular-momentum supported 
disk sizes. There is long history of studies relying on 
spin parameter distributions derived from dissipationless 
simulations to predict disk sizes and bulge properties in 
galax i es Ce.g..lFall fc Efstathioulll980l: iBlumenthal et ail 
MM IMo et all 119981 : IBullock et all 120011 ). Given that 
cool halo gas is the primary compone nt that actually 
builds galaxies (e.g., iDekel et"aLll2009f ). our result that 
A coo ig as ~ 4Adm motivates reconsideration of the is- 
sue of disk sizes in a cosmological context (see also 
IMaller fc Dekell2002l) . 

Before moving on, we note that while we have included 
disk gas and halo gas in our A determinations, the gas as- 
sociated with the central disk typically constitutes only 
~ 10 — 20% of gas within this radius, and is therefore 
unimportant to the quantities plotted here. That is, the 
spin parameters we are discussing are dominated by halo 
gas. We will present a more detailed study of the angu- 
lar momentum buildup of these galaxies in an upcoming 
paper. 

The relatively large amount of angular momentum in 
the cool component suggests that cold-accretion material 
may display a smoking-gun signature in its kinematics. 
In the next section we demonstrate not only that this gas 
is rapidly orbiting about the halo center, but that it tends 
to rotate with the central disk. This is not surprising 
given that this is the gas that fuels disk formation in the 
end. What is important for our purposes is that this co- 
rotation between cool halo gas and the galaxy disk can 
be observed via halo absorption studies. 



4. KINEMATICS OF COOL HALO GAS 
4.1. Co-rotation of Cool Accreted Halo Gas 

The kinematics of cool gas in the halos of galaxies serve 
as an important probe of galaxy formation. Unfortu- 
nately, the detailed nature of cool halo gas is difficult to 
determine empirically, as observations are presently lim- 
ited to studying absorption along quasar lines of sight. 
As a result, even when the resulting sight line is close 
enough to a foreground galaxy to probe its gaseous halo, 
these observations only allow for a single absorption line 
of sight per galaxy. Fundamental properties such as the 
covering fraction and kinematics of cool halo gas can 
only be accumulated in a statistical manner, for many 
different galaxies, with varied properties, redshifts, and 
distances between the absorbing gas and the galaxy. De- 
spite these inherent difficulties, observations of cool halo 
gas have progressed substantially in the past decade — 
typically utilizing MgTI absorption, which arises in low 
ionization, metal enriched gas of fairly high density and 
probes neutral hydrogen column densities of ~ 10 16 ~ 22 
cm' 2 (IRiebv et al.ll2002t IChurchill et alJl2000h . 

Still, it remains unclear what fundamental processes 
give rise to the absorbing gas. Different observa- 
tions suggest that i t ma y result from extended thick 
disks (jSteidel et al.1 12002D. cool clo ud formation from 
hot halo gas (jChen fc Tinker! 120081 ), cool gas that is 
outflowing from ga l axy-scale winds dBond et al.l [20011: 
Bouche et ail 120061: [Prochter et all 120061: iWeiner et al.1 
20091 : iSteidel et all 120101: iRubin eTaTT l2011h. tidal 



streams from recent m erger activity ( Churchill et al.l 
l2005tlRubin et aD2010| ). or a complex combination of the 
above, wit h inclusion of cosmolqgical g as accretion along 
filaments (|Kacprzak et al.l 120071120101) . In this section, 
we create mock observations of our simulated galaxies, 
utilizing line of sight absorption profiles of cool gas (as 
detailed in £12. 2j) to probe observable kinematic signatures 
of cool accreted gas. 

Before analyzing our mock absorption lines, we give 
a visual impression of the rotation of cool gas within 
R < 100 co- moving kpc in Figure [5] For each pixel in 
the image, the blue-red shading corresponds to the mass- 
weighted average velocity of HI gas along the line of sight. 
The galaxy, seen near edge-on in the inner region of the 
figure (at an inclination of 75 degrees), shows a clear 
rotation signature; gas on the left side of the disk shows 
a negative (blue) line of sight velocity while the right side 
of the disk has a positive (red) line of sight velocity. As 
shown, the rotation signature extends well beyond the 
galactic disk; in fact, the majority of cool halo gas out to 
R < 100 co-moving kpc is co-rotating with the galaxy. 

This kinematic correlation is not surprising, given that 
the angular momentum vectors of the inner and outer re- 
gions of dark matter halos are expected to be in rough 
alignment (typical misalignment ~ 30 degrees ; e.g. , 
van den Bosch et~aLl 127)021 ISh arma fc Steinmetzl 127)051 : 



Bailin et al.l l2005t iBett et all 120101 ). Nor is it surpris- 
ing that cold mode gas has high angular momentum 
before falling in to build the disk, as multiple past 
investigations into the properties of cold mode accre- 
tion have noted the presence of high angular momen- 
tum cool g as in proximity to (but more extended than) 
the galaxy (iKeres et al.ll2009bl: [Keres fc Hernauistl l200S 



lAgertz et ai.1 120091 : iBrook et al.H2011f ) Indeed, the grav 
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Fig. 6. — Co-rotation of cool halo gas in Halo 2's main galaxy at z — 1.4. The large black circles denote a radius of 100 co-moving kpc 
and the logarithmic color scale shows the projected column density of HI. In the left panel, observable sight lines (-/Vhi > 10 16 cm ) are 
categorized as orbiting (filled circles) or not orbiting (open squares). In the right panel, the direction of the velocity offset of orbiting sight 
lines is explicitly displayed, with redshifts and blucshifts given by circle-X and circle-dot symbols, respectively. 



itational potential energy of infalling gas should result 
in velocities of ~ 100 — 300 kms -1 , and if this energy 
is not dissipated by shock-heating, the conversion of po- 
tential energy into rotation al velocity is a qui te natural 
possibility (see discussion in iKeres et all [2001 . 

While Figure [5] gives an intuitive visual impression of 
this rotation, a mass weighted velocity of cool gas along 
a line of sight is not directly observable. In order to de- 
termine if this co-rotation signature might be observed 
via studies of absorption systems, we analyze 794 regu- 
larly spaced absorption profile sight lines within R < 100 
co-moving kpc of the central galaxjQ. Comparing each 
absorption profile to the rotation signature of the galactic 
disk, we categorize sight lines as orbiting or not orbiting 
based on their velocity offsets, and further divide orbiting 
sight lines into those that are co-rotating or anti-rotating 
with the galactic disk (see ^2.2l for detailed definitions of 
these classifications). 

Figure [6] again shows Halo 2 at z = 1.4 (this time 
viewed edge-on), but presenting a more observationally- 
oriented method than before. For both panels in Figure 
[6l the image width is 200 co-moving kpc, with the black 
circle marking a 100 co-moving kpc radius from the cen- 
ter of the galaxy. The colored contours in the background 
of each panel show the projected column density of HI 
with the color code on a log scale from > 10 16 cm" 2 
(red) to > 10 22 cm" 2 (purple/black). On this scale, the 
edge-on disk is visible as the dark horizontal bar across 
the very center of the image. We remove sight lines in 
the central region (R < 10 kpc) from our analysis, since 
we are interested in the properties of the cool halo gas, 
not gas that is likely associated with the galactic disk. 

13 The 794 sight lines form a square grid of 32 2 cells within a 
cube of length 200 co-moving kpc per side, after removing cells 
with R > 100 co-moving kpc. 



(This is why there are no symbols in the center of ei- 
ther image.) Along this projection, the central galaxy is 
rcdshifted on the left and blueshiftcd on the right. 

In the left panel, we have overlaid symbols representing 
the division between sight lines with significant velocity 
offsets (orbiting sight lines; filled circles) and those that 
do not (open squares). As the figure makes clear, the 
vast majority of sight lines are orbiting about the central 
galaxy, with velocities offset from the systemic velocity 
of the system. At this particular snapshot, the orbit- 
ing fraction (fraction of all iVni > 10 16 cm -2 sight lines 
which are orbiting) is ~ 90%. 

The right panel is identical to the left, except that 
orbiting sight lines (filled circles) now have their veloc- 
ity offset shown explicitly, so that we may determine 
which are co-rotating or anti-rotating with the galaxy. 
Sight lines that are redshifted with respect to systemic 
("into" the image) are shown by circle-X symbols, while 
blueshifted sight lines ( "out of" the image) are given by 
circle-dots. For this example snapshot and orientation, 
co- rotating sight lines follow the kinematics of the galaxy: 
redshifted on the left side of the image or blueshiftcd 
on the right side of the image. The co-rotation fraction 
(fraction of -/Vhi > 10 16 cm -2 sight lines that show this 
co-rotation signature) is ~ 70% for this imagcF^I. 

4.2. Effects of Galaxy Inclination 

In order to investigate the generality of our results, we 
repeat our kinematic analysis from §4.11 for four equally 
spaced viewing angles around each galaxy (at the same 
inclination, from — 135 degrees, inclusive) as well as 

14 Since lines of sight perpendicular to the plane of the galaxy 
cannot be classified as co-rotating or anti-rotating, we leave these 
as open circles, and do not include them in computation of the 
co-rotation fraction. 




Fig. 7. — Observable characteristics of halo absorption line kinematics as a function of central galaxy inclination angle for Halo 2 at 
z = 1.4. At each angle, we have analyzed four viewing orientations around the galaxy. Left and Middle: the orbiting and co-rotation 
fractions for different galaxy orientations. Squares and triangles show R < 50, 100 co-moving kpc statistics, with each symbol representing a 
different orientation of the galaxy. The orbiting fraction is universally high, while co-rotation fractions vary strongly with galaxy inclination. 
Bight: magnitude of velocity offset with galaxy inclination. Each curve shows the relative frequency of different velocity offsets for different 
galaxy inclinations, for lines of sight with R < 100 co-moving kpc and TVjji > 10 16 cm* 2 . The vertical lines indicate the median velocity 
offset for each inclination angle. 



seven equally spaced galaxy inclinations (0 — 90 degrees, 
inclusive) for a total of 28 unique orientations per galaxy. 
Continuing with the z = 1.4 snapshot of Halo 2, Figure 
[7] shows how the orbiting fraction (left panel) co-rotation 
fraction (middle panel) and velocity offset distributions 
(right panel) vary with these different orientations. In 
the left two panels, the x-axis shows the 7 different incli- 
nations of the galaxy, while each of the 4 symbols along 
the y-axis shows one of the 4 different viewing angles 
around the galax-yP^l. In order to explore how the orbit- 
ing and co-rotation fractions might vary with radius, we 
show statistics for two choices of outer radii (squares, tri- 
angles for R < 50, 100 co-moving kpc, respectively). The 
dashed and solid curves gives the mean values at each 
inclination for R < 50 and R < 100 co-moving kpc. 

Regardless of the galaxy's orientation, the overall or- 
biting signature stays constant, with a mean orbiting 
fraction of ~ 70 — 90%. This velocity offset signature is 
not a strong function of radius from the central galaxy. 
Cool halo gas at R = 100 co-moving kpc seems equally 
likely to be orbiting about the galaxy in some fashion as 
gas that is closer in to the galaxy. 

In contrast, the co-rotation fraction varies strongly 
with the orientation of the galaxy. At a fixed inclination, 
the co-rotation signature shows scatter of ~ 20 — 30%, 
depending on the viewing angle, but even more drastic 
is the dependence on inclination. Averaging over the 4 
viewing angles, ~ 70% of absorption sight lines appear to 
be co-rotating with the galaxy when viewed near edge-on. 
When viewed near-face on, not only does this co-rotation 
signature drop, but the cool halo gas actually appears to 
be anti-rotating. 

The right panel of this figure investigates this issue fur- 
ther, showing the distribution of velocity offsets for dif- 
ferent galaxy inclinations (averaged over the four viewing 
angles). When the galaxy is viewed near edge-on (solid 



black line) , the cool halo gas shows a wide distribution of 
velocities, with velocity offsets up to ~ 300 kms" 1 from 
the galaxy's systemic velocity, and an average velocity 
offset of ~ 100 kms -1 (the vertical line). While most 
sight lines still show velocity offsets at higher inclinations 
(i.e., most sight lines are still orbiting) the magnitude of 
these offsets are smaller, typically within ~ 100 kms -1 
from systemic, with an average offset of ~ 50 kms" 1 . 
These smaller magnitudes allow for projection effects to 
grow increasingly important for lines of sight to near face 
on galaxies. 

In order to fully understand this trend with inclina- 
tion, Figure [3] presents six different orientations of this 
galaxy (the same galaxy as Figures [SJ El and [7] Halo 
2 at z = 1.4). The contours and symbols in Figure [8] 
are identical to Figure [BJ with the left, middle, and right 
panels showing th e g alaxy at inclinations of 90, 45, de- 
grees, respectiverjQ. The top and bottom panels show 
two different viewing angles around the galaxy (90 de- 
grees apart). 

Apparent in the upper left panel of this figure is that 
most of the cool halo gas around this galaxy is confined 
to a single plane: a warped disk of cool accreted gas 
that extends to a radius of ~ 50 — 75 co-moving kpc 
(not an unusual occur rence hydrodynam ic simulations; 
e.g. lAgertz etHl l2009; Ros kar et alJl2010T) . These struc- 
tures, which we choose to call cold flow disks, contribute 
a significant portion of the total cool halo gas in our 
galaxies and tend to align roughly with the large-scale 
inflow filaments (see Figure [T]) . Cold Sow disks tend to 
be a common and prolonged phenomena in our simula- 
tions, and will arise naturally in systems with the high 
specific angular momentum content we measure. 

The top left panel demonstrates that the cold flow disk 
is, in fact, co-rotating with the galaxy, since the vast ma- 
jority of absorption sight lines show co-rotation (circle-X 



15 The order of rotation for all 28 orientations is first by angle, 
then by inclination. 



16 The face-on panels at the right show a "top-down" view of 
the edge-on panels at the left. 
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Fig. 8. — Identical to Figure [6] but showing the effect of galaxy inclination on the co-rotation signature. The left, middle, and right 
panels show galaxy inclinations of 90, 45, degrees, such that the right panels are top-down views of the left panels. The top and bottom 
panels represent different angles around the galaxy (at fixed inclination, 90 degrees apart). Note from the top left panel that most of the 
cool co-rotating halo gas is in an extended warped disk structure. Due to projection effects, gas that is co-rotating in this extended disk 
appears to be anti-rotating when the galaxy is near face-on. 



on the left, circle-dot on the right); i.e., the warped disk 
is rotating in roughly the same manner as the galaxy. 
Now consider the projection effects of viewing this same 
galaxy from the top-down. Because the left side of the 
gaseous disk is now projected out of the image, and the 
right side of the gaseous disk is projected into the image, 
the top half of the image is blucshifted, and the bottom 
half is redshiftcd. Even though the gaseous disk is co- 
rotating with the galaxy, the warp of the extended disk, 
in combination with projection effects result in a per- 
ceived co-rotation fraction of ~ 50% when the galaxy is 
near face-on. 

The bottom panels show an even more severe exam- 
ple of projection effects. In the left image, the warped 
disk is now seen in projection along the image, into the 
image along the top half and out of the image along the 
bottom half. Again, the co-rotation signature is evident, 
with the left side of the warped disk being redshifted and 
the right side blucshifted, the same as the galaxy. Once 
this image is viewed face-on (bottom right panel), the 
projected warp of the disk has completely flipped; the 
top half is now projecting out of the image, and bottom 
half is projecting into the image. As a result, the line 
of sight velocity of the gas changes sign as well, and the 
entire warped disk appears to anti-rotate with the galaxy. 

As evidenced by these examples, projected co-rotation 



of cool halo gas via absorption is not a reliable indicator 
of true co-rotation unless the galaxy is viewed near edge- 
on. Even if 100% of cool halo gas is in rough co-rotation 
with the galaxy's disk (e.g. in a cold flow disk), absorp- 
tion lines around moderately face-on galaxies may give 
false impressions of anti-rotation. 

5. EVOLUTION OF THE ORBITING AND CO-ROTATING 
SIGNATURES 

In SJU we focused on a single example: Halo 2 at 
z = 1.4. For this galaxy at this particular point in time, 
we found an overwhelming signature of orbiting gas (cool 
gas completely offset from the systemic velocity if viewed 
in absorption) and co-rotating gas (cool gas rotating in 
the same direction as the galactic disk when viewed in 
absorption) , using mock absorption profiles along lines of 
sight through the galaxy halo. Our goal has been to aid 
direct comparison to observations that utilize metal line 
absorption along lines of sight to backgrounds quasars. 
We now repeat our kinematic analysis for both simula- 
tions, and at a variety of redshifts, in order to deter- 
mine if our previous example was an unusual case, or 
if it points to a predictable signature in the behavior of 
cool accreted halo gas in galaxies. Because we found that 
projection effects result in unreliable rotation signatures 
unless the galaxy is viewed near edge-on (in regards to 
co-rotation versus anti-rotation; see ^4.2|) . we now limit 
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Fig. 9. — Covering fraction, orbiting fraction and co-rotation fraction versus time for both simulated galaxies, when viewed near edgc-on 
(inclination > 70 degrees). The solid and dashed curves are for Halo 1 and Halo 2, respectively. The error bars show the la scatter about 
the mean when averaging over 12 different viewing angles (see text). 



our analysis to inclinations of i > 70 degrees. 

Figure [5] shows the covering fraction (left panel), or- 
biting fraction (middle panel) and co-rotating fraction 
(right panel) of Halo 1 (solid red lines) and Halo 2 
(dashed black lines) for sight lines within R < 100 co- 
moving kpc, and with TVhi > 10 16 cm~ 2 . Error bars 
show the a scatter due to different galaxy orientations, 
as discussed in §4.21 In the right two panels, the curves 
change from thick to thin (and we discontinue error bars) 
when the galaxy transitions from cold mode accretion to 
hot mode accretion (sec £16.31 for more discussion of this 
transition). After this transition, the covering fraction 
of absorption gas drops precipitously (as seen in the left 
panel), resulting in small number statistics and large er- 
ror bars. We also note that the total mass in cool halo 
gas drops significantly after this transition as well (Fig- 
ure 2|) and our galaxies evolve into more bulge dominated 
systems (due to major mergers and secular processes, 
combined with the lack of a reservoir of cool high an- 
gular momentum halo gas) making the very definition 
of co-rotation with a presumed galactic disk nonsensical. 
Though not shown in this figure (for clarity), we find that 
the results in the two right panels are nearly identical for 
sight lines within R < 50 co-moving kpc (as opposed to 
R < 100), indicating a very weak dependence on radius, 
if any. 

Although our statistics are limited to a single pair of 
simulated galaxies, the middle panel of Figure |9] suggests 
that very high orbiting fractions (~ 80 — 95%) may be 
a natural consequence of cosmological gas accretion in 
LCDM, even for galaxies with quite different merger his- 
tories, such as our two galaxies. The co-rotation frac- 
tion, though, shows distinct differences between our two 
galaxies. Halo 2 experiences a quiescent merger history 
at z > 0.8, and shows a very strong signature of co- 
rotation once it builds a stable disk galaxy at its center 
(~ 60 - 80% of all sight lines for z < 2.5). Halo 1 has 
a similar accretion history (and co-rotation fraction evo- 
lution) until z 2.7, at which point it experiences a 
quick succession of two large mergers at z ~ 2.3. While 
mergers at high redshift are commonplace, the infall tra- 
jectory of the merger at z ~ 2.3 is strongly misaligned 



with the overall rotation of the galaxy (and halo) at that 
time. As a result, this merger completely re-orients the 
angular momentum axis of the galaxy, so that Halo l's 
galaxy at z = 1.5 is rotating in the opposite direction as 
its progenitor at z = 2.5. Without knowing how often 
major mergers result in such a drastic redistribution of 
angular momentum, it is impossible to speculate on how 
representative Halo 1 is, in terms of the kinematics of 
its cool halo gas between z ~ 1.5 — 2.5. Even in this se- 
vere case, with the angular momentum axis of the galaxy 
changing drastically due to a large retrograde merger, if 
we disregard a brief transition period near z = 2, Halo 1 
still shows an average co-rotation fraction of ~ 50 — 60% 
of all sight lines. 

6. OBSERVATIONAL IMPLICATIONS 
6.1. Accretion Signature Versus Outflows 

Summarizing our results from $5j wc find that ~ 80 — 
95% of absorption sight lines in our simulated galaxies 
orbit about the central galaxy, with significant velocity 
offsets from systemic. The co-rotation fraction of all sight 
lines varies in our two galaxies, based on their merger 
histories, from - 50 - 60% for Halo 1 to - 60 - 80% 
for Halo 2 (z < 2.5). Although our lack of radiation 
transfer may result in quantitative changes to these re- 
sults, the qualitative kinematic behavior of our galaxies 
cool gaseous halos appears independent of column den- 
sity. As such, the orbiting and co-rotating signatures 
we have presented here should be observable via halo 
absorption systems. We note, however, that since the 
co-rotation signature require galaxy kinematic and incli- 
nation information, this signature is less useful for high 
redshift observations, where such information is difficult 
(if not impossible) to acquire. 

Given the large amount of angular momentum asso- 
ciated with the accreted gas in our simulations, it is 
hard to imagine a scenario where radial/bi-conical out- 
flows from galactic disks would show such strong orbiting 
signatures. In principle, gas rotating in the disk could 
be blown directly upward. In the absence of angular- 
momentum robbing interactions with the ambient hot- 
halo gas, this gas would preserve a sense of rotation 
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in near-polar projections from the disk. However, the 
near-polar projection is precisely where infalling gas is 
least likely to show velocity offsets. One would expect 
these outflows to be detected along precisely those pro- 
jections where infalling gas is least likely to show velocity 
offsets — where the co-rotating signature is most affected 
by projection effects for a near edge-on disk-near perpen- 
dicular to the plane of the galaxy. Of course, perfectly 
polar sight-lines in such a scenario would also show two- 
sided velocity peaks (or a broadened peak about zero). 
This is again quite different from the signature expected 
from accreted gas. 

Moving away from the poles, the expected signal for 
rotating accreted gas is that a large fraction of the sight 
lines within ~ 45 degrees of the disk plane should show 
velocity offsets of ~ 80 kms" 1 or more at projected dis- 
tances of 100 co-moving kpc (40 physical kpc at z=1.5) 
from the galaxy (see Figure [5]). Consider gas that ini- 
tially samples the angular momentum in a disk rotating 
at ~ 200 kms^ 1 and with a half-mass radius of ~ 5 kpc. 
If that gas were blown out to projected radii of 40 kpc 
without losing any angular momentum to hot-gas halo 
interactions, we would expect typical velocities of only 
~ 25 kms , which is significantly smaller than what 
can be expected for the accreted gas. Furthermore, this 
estimate may be considered maximal in the sense that 
outflows are expected to preferentially expel low angular 
momentum material from galaxies, rather than ejecting 
disk gas that uniformly samples velocities within the disk 
(IMaller fc Dekelll2lM liTrook et alJl20Tl . 

One possible concern is that the inclusion of outflowing 
gas may radically disrupt the behavior of the accreted 
gas, such that it no longer shows a clear signature of 
orbiting/co-rotating about the galaxy. While it is be- 
yond the scope of this paper to investigate the effects 
of outflows on cold m ode gas accretion, a recent study 
by iBrook et al.l ()201lD found that while outflowing gas 
in their simulations (which tended to be bi-conical, near 
perpendicular to the disk) had a noticeable impact on 
the overall dynamics of their galaxies, the high angular 
momentum extended gas disks were least affected by the 
outflowing gas, since these extended disks arc (roughly) 
in the plane of the disk and outflowing material is blown 
perpendicular to the disk. 

6.2. Galaxy-Absorber Kinematics 

In order to compare our findings to observations, we 
combine results from several metal line absorption stud- 
ies th at compare the kinematics of the galaxy-absorber 
pairs dSteidel et alj|2003 lEllison et al.|[200l iChen et all 
120051: iKacprzak et al.ll2010f) . The redshift range for the 
total combined sample is 0.1 < z < 1.0, with typical 
projected radii R ~ 20 — 100 physical kpc between the 
galaxy and the absorber, and typical galaxy luminosi- 
ties of ~ 0.6L*. While our simulated galaxies grow too 
massive at z < 1 to make a direct comparison in the 
same redshift interval (their cold mode accretion ends, 
and they quickly lose their cool gaseous halos) our galax- 
ies are of comparable luminosity ~ 0.4L* at z > 1. While 
we urge the reader to keep this redshift discrepancy in 
mind, we believe a useful comparison can still be made, 
since the covering fraction of infalling cool gas appears to 
be more strongly corre lated with halo mas s than redshift 
(see Figure [9] also see lStewart et al.ll2011l) . 



Limiting these observations to galaxies with disk- 
like morphology (and not counting repeat analysis of 
the same object), these four studies contain (5,1,3,10) 
unique pairs of objects, respectively. Adopting our same 
definitions of orbiting and co-rotating sight lines as de- 
fined in §2.2[ we find that (4, 1, 2, 7) of these absorbers 
are orbiting, with an observed orbiting fraction of 74% 
(14/19). 

Of the 14 orbiting systems from this combined sample, 
(1, 0, 0, 1) systems were at a near right angle to the galac- 
tic disk. These ambiguous systems are removed from 
further comparison. Classifying the remaining 12 sys- 
tems, we find that (3,1,2,3) showed co-rotation, with 
only (0, 0, 0, 3) showing anti-rotation. This suggests that 
roughly 75% (9/12) of orbiting systems are co- rotating, 
for a total co-rotation fraction of 56%. A summary com- 
parison of these observed values to our simulations is 
shown in Table [TJ including Poisson \fN errors on the 
total number of galaxy-absorber pairs used to compute 
each statistic. In general, we find that these observations 
are in very good agreement with our predicted values, 
suggesting that these absorption studies may constitute 
the first indirect observations of cosmological gas accre- 
tion onto galaxies. Indeed, a comparison to cosmological 
simula tions was part of the analysis of IKacprzak et all 
(|2010D . They concluded that the absorption systems 
were consistent with the properties of cosmologically in- 
falling gas in LCDM. 

We also compare to IKacprzak et al.l (|2011| ) , who per- 
formed a similar analysis for a sample of absorber-galaxy 
pairs around more massive galaxies at z ~ 0.1 for more 
luminous galaxies ~ 1.9L*. Only 54% (7/13) of sight 
lines in this sample orbit the galaxicEFl. The lower or- 
biting fraction suggests that these absorption systems 
show more complex kin ematics than fres h ly acc reted gas, 
though the analysis of iKacprzak et al.l (|2011[ ) suggests 
that they are not well fit by purely outflowing gas, ei- 
ther. For purposes of co-rotation vs. anti-rotation, we 
remove from their sample 2 absorption lines near perpen- 
dicular to the slit across galaxy, as well as 2 sight lines 
around elliptical galaxies that are not dominated by or- 
dered rotation. Of this modified sample (which contains 
3 orbiting systems) only 1 sight line shows co-rotation 
with the galaxy. 

This co-rotation fraction is slightly lower than expected 
for galaxies above the critical transition mass for shock- 
heating infalling gas, though the qualitative trend agrees 
with our findings. Indeed, their sample consists of galax- 
ies with SDSS M r magnitudes in the range (-20 to -21), 
appropriate for those inhabiting dark matter halos with 
masses t hat span the transit ion mass M v ; r ~ 10 12 A/© 
(see e.g., iTollerud et all 1 20 111). Unlike the les s massive, 



higher redshift samples ( Kacprzak et al.ll2010l as well as 
Steidel, Ellison, Chen, and collaborat ors) the massive, 
low-z galaxies of IKacprzak et aLl (|2011l ) are not expected 
to have halos dominated by freshly-accreted (and orbit- 
ing) cold mode gas. The paucity of cold mode accretion 
for massive galaxies means that the cool gas present in 



17 IKac przak et al. ( 2011]) quotes a velocity offset fraction of 5/13, 
rather than 7/13, because two of their galaxies have most of their 
absorption offset to one side, but not all. We include these two 
cases as orbiting for the sake of this comparison, but neither defi- 
nition changes our qualitative conclusions here. 
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TABLE 1 

ORBITING AND CO-ROTATING SIGNATURES: 
SIMULATIONS VS. OBSERVATIONS. 





Observations^ 


Halo 1 


Halo 2 


Redshift range 


0.1 - 1.0 


1.3 - 


- 2.5 


0.J 


- 2.5 


Orbiting 


74 ± 20% 


80 - 


95% 


80 


-95% 


Co-rotating 


56 ± 18% 


50 - 


60% 


60 


-80% 



t Combined sam ple of Steidel ct al. (2002); Ell ison et all lf2003T ): 
IChen et al.l pOOlih ; IKacprzak et al.l l|2010l ) 

these systems is more likely to be the result of outflow 
material which should not display the co-rotation signa- 
tures outlined in this paper. Alternately, (unless it is 
associated visible infalling satellite galaxies) cool gas in 
these massive galaxies may also be the result of cloud 
condensation from the hot halo, which is roughly equally 
likely to be co-rotating or anti-rotating (see right panel 
of Figure EQ after the shock-hea ting transition). The 
analysis of lKacprzak et all (|2011D suggests that the sys- 
tems observed are not likely to be environmental galaxy- 
galaxy effects, nor from star formation driven winds, 
making gas accretion from the hot halo the most likely 
option for this sample. 



6.3. Covering Fractions 

So long as galaxies experience cold mode accretion, 
infalling gas results in moder ate covering fractions the 
inner halos of these galaxies (|Faucher-Giguere fc Keresl 
120111: iKimm et all 120111 : iFumagalli et all 120 111 ), though 
the resulting covering fractions are smaller than ob- 
served for high redshi ft galaxies domina ted by signa- 
tures of o utflows fe.g.. ISteidel et al.ll2010D . In a previ- 
ous work (jStewart et al.l 120111) . we examined the cov- 
ering fraction as a function of halo mass for the same 
two simulated galaxies studied in this paper, finding 
that once each galaxy crosses a critical threshold in 
halo mass, M v i r ~ 10 12 M Q , almost no cool accreted 
gas reaches the inner region of the halo. The under- 
lying cause for this transition is because beyond a cer- 
tain mass, halos are capable of sustaining shocks in in- 
falling gas, allowing the vast majority of cool accreted 
gas to shock- heat fsee e.g., [Birnboim fc Dekell 12005 
Keres et al.l 120051: iDekel fe Birnboiml l2006TlKcrcs et al. 
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Faucher-Giguere et al.l 12011: | van de Voort et al.l 
2011[j fe l. As a result. iStewart et al.l (|2011l) found a strong 
trend between the covering fraction of cool halo gas and 
halo mass; f c (< 50 co-moving kpc) ~ 30 — 50% before 
this transition, and only ~ 5% afterwards. This sharp 
decline in covering fraction is also shown in the left panel 
of Figure [9] 

This transition in covering fraction should be observ- 
able in absorption studies, provided there is some way 
to distinguish accreted gas from other sources of cool 
halo gas (e.g. outflows from the galaxy). The results we 
have presented here allow observations to do precisely 
this. For any given absorption line through a galaxy's 
gaseous halo, cosmologically accreted cool gas should al- 
most always orbit about the galaxy (orbiting fraction of 

18 The critical halo mass where galaxies begin to accrete more 
hot gas than cool gas is considerably lower. Using this definition, 
the transition mass is often quoted as ~ 10 11,5 Mq, rather than 
1O 12 M . 



80—95%), and most sight lines should be co-rotating with 
the nearest side of the galactic disk, as long as the galaxy 
is viewed near edge-on (co-rotation fraction ~ 60—70%). 

Coupling the res ults of this paper with those of 
IStewart et al.l (j 2 1 If ) . we predict a clear, observable sig- 
nature in the kinematics of galaxy-absorber pairs as a 
function of halo mass. Below the critical halo mass, some 
portion of the cool halo gas detected in absorption will 
be due to cool accreted gaj 1 ^!. This accreted gas will co- 
rotate with the galaxy. 

Above the critical halo mass, cool accreted gas no 
longer contributes substantially to the covering fraction, 
so even if the total covering fraction remains high (as a 
result of outflows) the absorbing gas will no longer co- 
rotate. As a result, the covering fraction of co-rotating 
sight lines should drop dramatically at or near this critical 
halo mass. While current observations lack the statistical 
power and halo mass range to robustly test this predic- 
tion, we believe that future observations that continue to 
compare the properties of many absorption systems with 
their associated galaxies will be vital in furthering our 
understanding of gas accretion onto galaxies, as well as 
the nature of gaseous halos around galaxies. 

6.4. Connections with Extended XUV and HI Disks 

The extended cold flow disk structures visual- 
ized in Figure El may help expl ain the existence 



HI 



of XU V disks (iThilker et alJ I2ll05l l200l. large 
disks (iGarcia-Ruiz et al.1 120021 : lOosterloo et al.l 120071 : 
I Walter et al.l 120081 ). and some stellar disks (jBarthj |2003) 
that extend to ~ 50 kpc (or more in some cases). It is 
generally hard to understand the existence of disks with 
radial extents as large as ~ 50 kpc in models that link 
disk galaxy sizes t o dark matter s pins of Arim ~ 0-04 
(iFall fc Efstathioul H9801 IMo et all 119981 : IBullock et all 
120011: iDutton fc van den Bosch! I2009f ). In these simple 
models, disk material is expected to extend to radii 
i?disk — Ai? v ir = 10 kpc for galaxy-size halos of virial 
radius i? v i r ~ 250 kpc. More extended gas disks are eas- 
ier to understand when they arise from cold-mode accre- 
tion with spin parameters as high as A ~ 0.2 (Figure |4|). 
Therefore, we expect XUV disks and extended HI disk 
configurations to be more common in halos less massive 
than ~ 10 12 Mq (or with maximum circular velocities 
smaller than ~ 165 km/s), though a detailed comparison 
between the cold flow disks simulated here and observed 
XUV and extended HI disks is beyond the scope of this 
work. 

6.5. High-redshift Limitations 

While the high orbiting fractions predicted here for in- 
falling cold mode gas only requires a systemic velocity for 
the galaxy, the co-rotation signature has a notable limita- 
tion in that it requires detailed galaxy information (rota- 
tion curve, morphology, inclination) in order to compare 
galactic rotation with the associated absorption system. 
While this is possible for low redshift galaxies (e.g., the 
combined sample discussed previously), it is extremely 
difficult at higher redshifts, making the co-rotation sig- 
nature increasingly difficult to observe at high-z. 

19 We note that without knowing the full effects of outflows, it 
is impossible to say how large a portion of the detectable cool gas 
halo will come directly from fresh accretion. 
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In addition, because we have focused on cool gas in 
galaxy halos, rather than filamentary gas that is not as- 
sociated with a particular galaxy/halo, the radial extent 
of this orbiting gas scales roughly with the galaxy virial 
radius (this is why we have focused on co-moving coor- 
dinates in our analysis). As a result, the radial extent 
of this infalling gas (found here to be ~ 100 co-moving 
kpc) corresponds to relatively small physical separations 
at high-z (~ 50 and ~ 33 kpc at z = 1 and z = 2). 
Once again, this is an added difficulty in observing this 
signature at high redshift, as it is difficult to find bright 
background objects within such small projected separa- 
tions to foreground galaxies. 

7. CONCLUSION 

We have used two high-resolution cosmological SPH 
simulations (with star formation and feedback prescrip- 
tions that have proven effective in reproducing realistic 
galaxies at z = 0) tool for studying the detailed 

properties of cool gas in galaxy halos. By creating ab- 
sorption sight lines probing the column density of cool 
gas in our simulated galaxies, we make predictions that 
should be comparable to observations, either via metal 
line absorption or direct detection of HI, though we note 
that without radiative transfer, the precise values pre- 
sented here are less robust than our qualitative findings. 
Our primary results are summarized as follows: 

1. Cosmologically accreted cold- mode gas contains 
significant angular momentum compared to the 
dark matter, and orbits about the galactic halo be- 
fore falling in to build the central disk. The dimen- 
sionless spin parameter of cool halo gas is ~ 3 — 5 
times higher than that of the dark matter, with 
A CO oigas ~ 0.1 — 0.2. Cold flows clearly supply not 
just baryonic mass to galaxies, but angular momen- 
tum as well. 

2. This orbiting cool halo gas should be observable via 
background object absorption line studies as lines 
that are offset from the galaxy's systemic velocity 
by - 50 - 100 km s" 1 . Specifically, - 80 - 95% 
of sight lines that show absorption (TVhi > 10 16 
cm~ 2 ) arc offset from the central galaxy's systemic 
velocity in a single direction. This high fraction 
of lines that are "orbiting" seems an unavoidable 
consequence of cosmological gas accretion onto disk 
galaxies. It is universal for both simulations, at all 
epochs, and is independent of the orientation of the 
galaxy. The orbiting gas consists of accreted cool 
gas from the cosmic web and subsequently falls in 
to help build the disk. 

3. Much of the orbiting cool halo gas commonly takes 
the form of a fairly thick, warped, extended disk in 
both our galaxies. As a result, the fraction of all 
sight lines that co-rotate with the galactic disk is 
also relatively high (~ 55% for Halo 1 and ~ 70% 
for Halo 2). This is a unique, observable signature 
of cosmological gas accretion, since there is no rea- 
son to expect outflowing gas to co-rotate in this 
fashion. We warn that due to line of sight projec- 
tion effects of the warped disk, this co-rotation sig- 
nature can only be reliably measured if the galaxy 



is near edge-on (i > 70 degrees). Less inclined sys- 
tems may produce false anti-rotation signatures. 

4. We compare our orbiting and co-rotating signa- 
tures to a combined observational sample of galaxy- 
absorber pairs for galaxies of similar luminosities 
(ISteidel et aLll200l lEllison et al.| [200l [Chen et al.l 
120051 ; iKacprzak et all 12010 ). finding rough agree- 
ment between our predictions and these observa- 
tions (Table [1]). These studies may provide some 
of the first (indirect) observational evidence for cos- 
mological accretion of cool gas onto galaxy halos. 

5. The covering fraction of cosmologically accreted 
cool gas within a given radius, f c (<R), is a strong 
function of R. At small radii, the covering frac- 
tion is near unity, decreases to negligible values 
at R > 100 co-moving kpc, and is well-fit by a 
power law. The slope of this power law depends 
on the minimum column density of absorbing gas: 
fc(Nm,< R) = {R/R )~ i3{Nhi \ with our galaxies 
well fit by i?o ~ 10 co-moving kpc and /?(> 10 16 
cm~ 2 ) ~ 0.6 — 0.8 (see Equations 1-3). 

6. The extended, warped, and thick cold-mode disks 
illustrated in Figure 8 may be associated with the 
extended UV and HI disks observed for galax- 
ies in the nearby un i verse dThilker et all 120051 
[20071 : IQosterloo et~aTl [20071: IWalter et alfl2003) . 
If so, we expect extended disk configurations to 
be more common for galaxies in halos less mas- 
sive than ~ 10 12 A/q (or with maximum circu- 
lar velocities smaller than ~ 165 km/s). We 
note that HI and XUV disks as large as those 
sometimes seen (~ 50 kpc) are quite difficult 
to understand if the available angular momen- 
tum for the gas mirrors that in the dark matter, 
as is often assumed in si mple LCDM- based mod- 
els of disk formation (e.g. | Fall &; Efstathioul 1 1980 : 
iBullock et al.ll200lUDutton fc van den Boschll2009l 
where i?disk — Ai? V i r and A ~ 0.04). If these disks 
are instead supported by cold accreted material 
(with A ^ 0.1 — 0.2) their extent and structure is 
more naturally explained. 

Our expectation is that cold mode accretion will re- 
sult in circum-galactic gas that orbits with fairly high 
angular momentum. Spherical outflows will be charac- 
teristically radial with low angular momentum, in stark 
contrast to infalling gas. While bi-conical outflows may 
give a velocity offsets under certain projections, similar 
to the accretion signature discussed here, the resulting 
outflows will exhibit primarily polar configurations; this 
is precisely where infalling gas is less likely to show veloc- 
ity offsets and the orbiting/co-rotating signature is least 
reliable. Given a sufficiently large sample of galaxy ab- 
sorber pairs, it should be possible to determine the frac- 
tion of orbiting and non-orbiting halo gas as a function of 
host halo mass, and to test the expectation that less mas- 
sive halos (below the critical shock mass) should prefer- 
entially host a higher fraction of cold-mode/orbiting gas 
than their more massive counterparts. 
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